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Abstract 

In this paper we describe the structure of extremal stationary spherically symmetric black hole 
solutions in the STU model of -D = 4, = 2 supergravity in terms of four-qubit systems. Our 
analysis extends the results of previous investigations based on three qubits. The basic idea facili- 
tating this four-qubit interpretation is the fact that stationary solutions in D = 4 supergravity can 
be described by dimensional reduction along the time direction. In this D = 2> picture the global 
symmetry group S'L(2,M)^^ of the model is extended by the Ehlers S'L(2,M) accounting for the 
fourth qubit. We introduce a four qubit state depending on the charges (electric, magnetic and 
NUT) the moduli and the warp factor. We relate the entanglement properties of this state to differ- 
ent classes of black hole solutions in the STU model. In the terminology of four qubit entanglement 
extremal black hole solutions correspond to nilpotent, and nonextremal ones to semisimple states. 
In arriving at this entanglement based scenario the role of the four algebraically independent four 
qubit SL{2,C) invariants is emphasized. 

PACS numbers: 11.25.Mj, 03.65.Ud, 03.67.Mn, 04.70.Dy 
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I. INTRODUCTION 



Recently striking multiple relations have been discovered between two seemingly unre- 
lated fields: Quantum Information Theory (QIT) and the physics of black hole solutions 
in String Theory-""^. Although the physical basis for this black hole qubit correspondence 
(or black hole analogy) is still to be clarified, it has repeatedly proved to be useful for ob- 
taining additional insight into one of the two fields by exploiting methods and techniques of 
the other-"—. The main correspondence founcii'^i^^-'i^ is between the macroscopic entropy 
formulas obtained for certain black hole solutions in supergravity theories and multiqubit 
and qutrit entanglement measures used in Quantum Information Theory. The basic rea- 
son for this correspondence is the occurrence of similar groups of symmetry in these very 
different contexts. On the stringy black hole side the groups in question are the global 
symmetry groups oi D = 4 classical supergravities, and on the QIT one the groups of local 
transformations for entangled subsystems not changing their multipartite entanglement. As 
far as physics is concerned an attempt has been made to understand these mathematical 
coincidences in terms of wrapped brane configurations giving rise to qubits-. 

Apart from understanding black hole entropy in quantum information theoretic terms 
the desire for an entanglement based understanding for issues of dynamics also arose. In 
particular in the special case of the STU model^^ it has been realizect- that it is possible 
to rephrase the attractor mechanism^^ as a distillation procedure of entangled "states" of 
very special kind on the event horizon. Such "states" for D = 4 extremal static spherical 
symmetric solutions are arising from more general ones of the foru^^ 

|^(r)) = {V®V® V)(53(r) ® S2(r) ® Si{r))\j), t = -. (1) 
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where r is the radial distance from the event horizon, Zj{T) = Xj(r) — iyj{T), j = 1, 2, 3 are 
the scalar fields, qi, and , I = 0,1,2,3 are the electric and magnetic charges occurring 



in the STU model^^. As we see these quantities are organized into a complex three-qubit 
state. This instructive notation clearly expresses the triality symmetry of the STU modeU^. 
Moreover, the classical symmetry group of the model (i.e. SL(2, M)^^) is manifested in this 
formalism by the fact that apart from the unitary matrices V, is lying on the SL{2, M)^^ 
orbit of the "charge-state" I7). The unitaries V®^ provide an embedding of the SL{2,'R)^^ 
symmetry group of this N = 2 supergravity model into GL{2, C)^^. 
The state of Eq. ([1]) has a number of remarkable properties^i^ii^. 

1. The three-tanglei^ T123, the unique triality and SL(2,C)^^ invariant 3-qubit entangle- 
ment measure based on Cayley's hyperdeterminantr^ii^, for is related to the macroscopic 
black hole entropy in the STU model as 

S = 7rv^ri23(|^)) = vrv/ri23(|7))- (5) 

2. The norm of with respect to the usual scalar product in with complex conju- 
gation in the first factor is the Black Hole Potential^ Vbh- 

3. The flat covariant derivatives with respect to the Kahler connection are acting on {tp) 
as bit flip errors on the qubits. 

4. For iJPS'-solutions and for non-BPS solutions with vanishing central charge^ \iIj{oo)) is 
a GHZ-state^^. For non-BPS solutions with non-vanishing central charge the corresponding 
states are graph-states known from QIT— . In this respect moduli stabilization is related to 
a distillation procedure of states with special entanglement properties at the event horizon. 

5. On the horizon bit flip errors on lip) are supressed for BPS solutions and for non-BPS 
ones they are not. The non-BPS solutions can be characterized by the number and types of 
bit-flip errors. 

6. After solving the equations of motion one obtains the attractor flow Zj{T) in moduli 
space. There is a flow |^/'(t)) associated to this one. For the non-BPS seed solution^^ it is 
possible to study how the distillation procedure unfolds itself— with the following result. In 
the asymptotically flat region we are starting with a \ip{0)) having 7 nonequal nonvanishing 
amplitudes and finally at the horizon we get a graph state \iIj{oo)) with merely 4 nonvanishing 
ones with equal magnitudes. 

7. The magnitude of the nonvanishing amplitudes of such "attractor states" is propor- 
tional to the black hole entropy. The relative phases of the amplitudes reflect the structure 
of the fake superpotential. 



3 



8. If we are starting with the very special values for the moduli corresponding to flat 
directions'^ this uniform structure at the horizon deteriorates^-, with the interpretation of 
errors of more general types acting on the qubits of the relevant attractor states. 

In addition to these interesting results based on three-qubit states there are ones which 
strongly hint at the possibility that for a complete understanding of STU black holes we 
have to embed our three-qubit states into four-qubit onea^. In particular one can generalize 
Eq.([T]) by also including the warp factor U{t) occurring in the static, spherically symmetric 
ansatz for the 4D space-time metric 

into a new state \x) defined a.s^ 

lx(r)) = e^W|V^(r)). (7) 

For the non-BPS seed solution^^ it has been shown that the 7 nonvanishing r dependent 
amplitudes of this state depending on the charges, the moduli and the warp factor satisfy 
a system of first order differential equations. This finding conforms with recent work done 
within the framework of the first order formalism for non-BPS solutions based on the so 
called fake superpotentia^. 

Moreover, within the realm of the more general class of stationary solutions it is well- 
known that the warp factor taken together with the NUT potential^^ a forms another 
S'L(2,M) doublet , a doublet with respect to the Ehlers group^^. Hence it is natural to 
suspect that for stationary solutions objects like |x(t)) are really four-qubit states in dis- 
guised form with the Ehlers group acting on a hidden extra qubit. The properties of these 
hypothetical 4-qubit states should account for the first order formalism hiding behind the 
integrability of the non-BPS fiow equations. 

Recent investigations clearly demonstrated that this should indeed be the case^l"— 
The key observation is that stationary solutions in D = A supergravity can be elegantly 
described by dimensional reduction along the time direction^-. In this picture stationary 
solutions can be identified as solutions to a D = 3 non-linear sigma model with target 
space being a symmetric space G/H with H non-compact. The property that is of basic 
significance for us is that the group G in this case extends the global symmetry group G4 
of D = 4 supergravity, by also incorporating the Ehlers S'L(2,M). In our specific case the 
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N = 2 STU model can be regarded as a consistent truncation of maximal N = 8, D = 4 
supergravity with G4 = -£^7(7), truncating to SL{2,'R)^^. Timelike reduction in the general 
case then yields the coset -E'8(8)/S'0*(16), or in the case of the STU truncation the one 
Ms = 50(4,4)/5L(2,M)^l We then expect the four copies of 5L(2,M)s giving rise to the 
group of local operations acting on the four qubits. Here the fourth qubit which accounts 
for the Ehlers group will then play a special role. 

The manifold A^s is the target space of the aforementioned sigma model. It has been 
proved^^ that for such symmetric target spaces stationary spherically symmetric black hole 
solutions can be obtained as geodesic curves on this pseudo Riemannian target space. Such 
geodesic curves are classified in terms of the Noether charges of the solutions. In the case of 
the STU model the coset representative V of our target space Ais and the related Noether 
charge can be written in the form reminiscent of a 4-qubit state^i^i^^. Moreover the line 
element on A^3 can be written in the fornt^I Tr(P^) which turns out to be just the quadratic 
4-qubit invariant, one of the four algebraically independent invariants characterizing 4-qubit 
systems^. It has also been observed^i^^ that the entanglement properties of such 4-qubit- 
like states seem to be related to the fact whether the extremal solution in question is BPS, 
non-BPS or non-BPS with vanishing central charge. Based on this finding the authors of 
this paper— mention that there might be a connection with issues concerning the black hole 
qubit correspondence, though in this field the D = 3 reformulation has never been used. 
(See Egs. (5.51)-(5.52) of that paper.) 

The aim of the present paper is to show, that using the D = 3 picture such 4-qubit 
interpretation indeed emerges naturally. Moreover, after establishing the desired connection 
we see that in this framework many aspects of the usual three-qubit interpretation can be 
understood in a nice and unified way. 

The organization of this paper is as follows. In Section II. we present the background ma- 
terial on the STU model and the basics of the D = 3 picture emerging after reduction along 
the time direction. In Section III. in a four-qubit notation we reconsider the usual Iwasawa 
parametrization of the physical patch of the pseudo- Riemannian manifold Ai^. This formal- 
ism is exploited in Section IV. where we describe the line element on A^3 as the canonical 
quadratic 4-qubit SL{2,C)^^ invariant. Here after a sequence of 4-qubit transformations 
(Hadamard gates, phase gates, and permutations) a very convenient realization for the "vier- 
bein" V is obtained. These transformations correspond to a special choice of basis in TcA^s 
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similar to the ones used in Ref .— rendering the " quaternionic vierbein" covariantly constant 
with respect to the spin connection. In Section V. we discuss the structure of conserved 
charges in an entanglement based framework. Here we see how our remarkable three-qubit 
state of Eq. ([T]) originates from the geometric data on A^3. As an important generalization 
we write down a generalization of Eq.([T]) for stationary solutions when the NUT charge is 
not zero. Section VI. is devoted to an analysis of the static, spherically symmetric solutions. 
Our treatment is based on the algebraically independent 4-qubit SL{2,<C)^'^ invariants. It 
is shown that in the language of QIT extremal solutions correspond to nilpotent, and nonex- 
tremal ones to semisimple 4-qubit states. Nilpotent states are the ones for which all of the 
four algebraically independent invariants vanish. This picture is dual to the usual charac- 
terization in terms of nilpotent orbits. Next in this entanglement based approach a study of 
the usual BPS and non-BPS solutions with vanishing central charge, and the non-BPS seed 
solution is given. These investigations culminate in establishing an explicit connection be- 
tween the results of Ref.~ and some standard ones on four-qubit entangled systems in QIT. 
Finally we present our conclusions and comments in Section VII. In an Appendix for the 
convenience of the reader we also included some background material concerning four-qubit 
systems. 

II. THE STU MODEL 

In the following we consider ungauged N = 2 supergravity in d = 4 coupled to n vector 
multiplets. The n = 3 case corresponds to the STU model. The bosonic part of the action 
(without hypermultiplets) is 



Here J-"^, and / = 0, 1,2. . .n are two-forms associated to the field strengths of 
n + 1 U{1) gauge-fields and their duals. 

The z"^ i = l,2...n are complex scalar (moduli) fields that can be regarded as local 
coordinates on a projective special Kahler manifold. This manifold for the STU model is 




(8) 
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[5'L(2,M)/f/(l)]^^. In the following we will denote the three complex scalar fields as 

= - iy\ J = 1, 2, 3, > 0. (9) 

With these definitions the metric and the connection on the scalar manifold are 

= rj, = Zi. (10) 



The metric above can be derived from the Kahler potential 

ir=-log(8yil/2Z/3) (11) 

as = didjK. For the STU model the scalar dependent vector couplings ReA/}j and 
ImAfij take the following form 

/ 2X1X2X3 — X2X3 —X1X3 —XiX2\ 

—X2X3 X3 X2 
— X1X3 X3 Xi 

\ -X1X2 X2 Xi 



vij = ReAfij 



(12) 



/i/j = ImAfjj = -ymvs 
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(13) 



Xl |2;i| X1X2 X1X3 

X2 X1X2 \z2f' X2X3 
|2 



2/11/21/3 

\X3 X1X3 X2X3 \Z3\"/ 

We note that these vector couplings can be derived from the holomorphic prepotential 



(14) 



Fix) 



x^ = (x°,xV), 



(15) 



via the standard procedure characterizing special Kahler geometry^. 

Our aim is to describe stationary solutions of the Euler-Lagrange equations arising from 
the Lagrangian of the STU model in a four-qubit entanglement based language. It is well- 
known that the most general ansatz for stationary solutions in four dimensions is 



ds' = -e^^{dt + uf + e-''''habdx''dx\ 



2U] 



(16) 



= dA^ = d{i\dt + uj) + A^), (17) 

where a, 6 = 1,2,3 correspond to the spacial directions. The quantities U , A^, uja and 
hab are regarded as 3D fields, i.e. the ansatz above corresponds to dimensional reduction 
io D = ?) along the timelike direction. In achieving this we have chosen the gauge such 
that the Lie-derivative of A^ with respect to the timelike Killing vector vanishes, and have 
chosen coordinates such that the isometry corresponding to this Killing vector is just a 
(time) translation. In this case the quantities in Eqs. fll6tiT7|) are merely depending on x", 
a = 1,2, 3. The ansatz for the gauge fields A^ refiects its decomposition to terms parallel 
{^^)-, and orthogonal {A^) components with respect to the timelike Killing vector™. 

After performing the dimensional reduction to D = 3 our starting Lagrangian of Eq. ([8]) 
takes the following form^"^ 

£ = £1 + £2 + £3, (18) 

where 

Ci = ~VhR[h] +dU A *dU + ^e-^^{da + f/d^^ - ^^rff/) A *{da + i^jd^-^ - i'' dlj), (19) 

£2 = Gqdz' A (20) 
£3 = ^e-2^/x,jrfe' A *de + \e-'''' [i'\dli - viKde) A K^O - ^JLd^)- (21) 
Here the new (axionic) scalars a and are coming from dualizing uj and A^ by2^ 

dli = uijdi^ - e^^fiu * {dA-^ + i^du) (22) 

da = e^^ * doj + iUli - lidi^. (23) 

Note also that here the exterior derivative is understood on the (generally curved) spatial 
slice with local coordinates x-' , j = 1, 2, 3. 

The dimensionally reduced Lagrangian £ can be written in the nice form of 3D gravity 
coupled to a nonlinear sigma model defined on the spatial slice with target manifold^- 
M3 = 50(4,4)/5L(2,M)^^ with the Lagrangian 

£ = -^VhR[h] + (7^„9„$™a'^$" (24) 

where C"™, m = 1, 2, ... 16 refers to the scalar fields: U, a, with 1 = 0,1,2,3 and 

j = 1, 2, 3. Here the line element on Ai-^ defines gmn as ds'j^ = gmn^"^^"' with the explicit 
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form 

i(isXi3 = Gij{z,z)dz'dz^ + dU^ + ^e-^^ida + iid^' - i'dlif 

Ludi'di' + ^l'\dlI - viKde){dlj - vjLd^)] ■ (25) 



In this paper we are only discussing the special case of stationary, weakly extremal so- 
lutions i.e. solutions when the spacial slices are flat^i^^. Single centered black holes with 
spherical symmetry are of this type. In this case the dynamics of the moduli are decou- 
pled from the 3D gravity and the metric ansatz can be chosen to be the form 



ds"" = -e'^idt + u) + e-^^{dr^ + r\de^ + sin^edip), (26) 

with the warp factor depending merely on r. Now the equations of motion are equivalent to 
light-like geodesic motion on M.-^ with the affine parameter t = K Since M.^ is a symmetric 
space there is a number of conserved Noether charges associated with this geodesic motion. 
The most important ones are the electric and magnetic charges and qi and the NUT 
charge A)^^'^'^. Static solutions are characterized by the vanishing of the NUT charge i.e. 
/c = 0. In this case the dynamics is described by the Lagrangian of a fiducial particle in a 
"black- hole potential" Vbh 

(dTJ\ ^ dy^ 
^) +G~^ + e'''VBHiz,z,p,q), (27) 

with the constraint 

Here the black hole potential Vbh is depending on the moduli as well on the charges. Its 
explicit form is given by 

An alternative expression for Vbh can be given in terms of the central charge of = 2 
supergravity, i.e. the charge of the graviphoton. 

Vbh = ZZ + G'^D.Z) (DjZ) (30) 

where for the STU model 
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Z = e^l'^W = e^/^(go + ^igi + Z2q2 + zsQs + ziZ2Z3P^ - Z2Z3p^ - ziz^p'^ - ziZ2p^), (31) 
and Da is the Kahler covariant derivative 

Da = {d, + ^m)z, (32) 

and W is the superpotential. 

Extremization of the effective Lagrangian Eq.f l27p with respect to the warp factor and 
the scalar fields yields the Euler-Lagrange equations 

U = e^^VsH, z' + r;.,i^i^ = e'^^d'VBH. (33) 

In these equations the dots denote derivatives with respect to r = ^. These radial evolu- 
tion equations taken together with the constraint Eq.( l28ll determine the structure of static, 
spherically symmetric, extremal black hole solutions in the STU model. For the more gen- 
eral stationary case with nonvanishing NUT charge the motion along , C,i and a does not 
separate from the one on U and z^ . In this case we obtain the generalization of Eqs. (l33|l . 
Since for our four-qubit picture we will not consider solutions of such kind we will not give 
the corresponding equations here. 

As we have seen from this section the radial evolution associated to stationary spherical 
symmetric black hole solutions of the D = 4 STU model can be described^"^ as geodesic 
motion in the moduli space A^3 of a dimensionally reduced D = 3 theory. The key issue of 
this reduction relevant to this paper is the enlargement of the D = 4 symmetry group from 
SL{2,Ry^ to the D = 3 one 50(4,4) containing SL(2,M)^^ as a subgroup. This result 
paves the way for the possibility to reinterpret our STU black holes as four-qubit systems. 

III. THE IWASAWA PARAMETRIZATION AND FOUR QUBITS 

Our starting point is the Iwasawa parametrization of the coset Ais = 5*0(4, 4) / 5*0(2, 2) x 
SO(2,2) ~ SO(4,4)/5L(2,M))^4 as used in the paper of Bossard et.al.-2I For this 
parametrization the 16 dimensional coset is (locally) coordinatized by the fields Xj,yj, 
4> = 2U,a, and the potentials and ^/ quantities featuring the Lagrangian £ of Eq.( l24|) . 
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In order to avoid using disturbing factors^ of a/2 we rescale the potentials and define new 
quantities , (i as 

= V2e, Ci = (34) 



In terms of these quantities the coset representative is 



(35) 



Here the four copies of SL{2, M) generators Ha, E^, F^, a = 0, 1, 2, 3 satisfy the commutation 
relations 

[Ea,Fa] = Ha, [Ha, Ea] = 2Ea, [Ha, Fa] = — 2Fq,, (36) 

and the 16 generators of so(4, 4) not belonging to the sl{2) © s/(2) © s/(2) © s/(2) algebra are 
denoted by the symbols Epi , Eqj, Fpi , Fqj, I = 0,1,2,3. This decomposition of generators 
answers the split 

so(4, 4) = [sl{2, M)]^ © (2, 2, 2, 2) = /i © m, (37) 

which we would like to explicitly describe. (For an explicit connection between our conven- 
tions described below, and the one as given by Bossard et.al.^^ we refer the reader to the 
Appendix.) 

The Lie-algebra so(4, 4) adapted to our 4-qubit description will be regarded as the set of 
8x8 matrices V satisfying 



where 



G 
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9. 



VG + GV^ = 



g = e ®e 



1 
-1 



T^{s3, S2, Si, So; D) 



An element of so(4, 4) will be parametrized as 

S3 © /2 + /s ® S2 Dg 

-D^g si © Jo + /i © So , 
Here the m-type generators are labelled by a real 4x4 matrix 

/ -Doooo -Doooi -Dooio -Dqoii^ 

-Doioo -Doioi -Doiio -Doiii 

-Diooo -Diooi -Dioio -Dioii 

\-Diioo -Diioi -Diiio -Dull/ 



D 



(38) 



(39) 



(40) 



(41) 
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which is expressed in terms of the amphtudes of a 4-qubit state with index structure 

Di3i2hio, i3,i2,ii,io = 0,l. (42) 

Notice that for convenience we have labelled the qubits from the right to the left. Moreover, 
the first qubit will be regarded as special explaining the somewhat unusual label: tQ. 
The h type generators are featuring the 2x2 matrices Sa of the form 

s^=\" , a = 0,l,2,3. (43) 

\fa J 

These matrices are expanded in terms of the ones 

-(::)• -(::)• -(::)• 

satisfying the relations of Eq.f l36|) . 

The labels of the 2x2 matrices appearing in Eq. ( |40|) are referring to the qubits they 
act on. This action is induced by commutators of the form [h, m] C m. More precisely after 
commuting the block off-diagonal m part with the block-diagonal h one using 

s£ + £s'^ = 0, sesl{2) (45) 

we get the action 

(s3 ^h + h® S2)D + D{sl ® Jo + /i ® si), (46) 
which is the first order term in the SL{2,M)^^ group action 

{S3®S2)D{Si®Sof, ^„G5L(2,M), a = 0,1, 2, 3. (47) 

Clearly this action in 4-qubit notation reads as 

Disi2hiQ^ 5Z ('S'3)i3i^,('S'2)j2i'2('S'l)ni',(5'o)ioiJ,A^,ii.i'ijJ,, (48) 

or in the notation used in Quantum Information Theory 

\D) ^ {S3 ® ^2 ® ^1 ® So)\D), \D) = J2 Di,^2^lio\Wlio)■ (49) 

*3*2*1*0=0,1 

We remark that for the convenience of the reader in the Appendix we included more details 
on the correpondence between the structure of the group 50(4, 4) and 4-qubit entanglement. 
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Now returning to our coset representative of Eg. (1350 . we introduce the new coordinates 



2U 



xo = cr, I/O = = e 



Using our 4-qubit realization in these coordinates we have 



(50) 



n 

0=0 



g-i logyaHa^-XaEa 



M3 (g) M2 

Ml (g) Mo 



(51) 



where 



M„ 



1 /I -Xo 



As a next step we introduce the 4x4 matrix and its associated 4-qubit state 



(52) 



c 



f -Co Ci o\ 
C2 

Ca 



(53) 



'^Coooo Coool Coolo Cooii\ 
Coloo Coioi Coiio Com 

ClOOO ClOOl ClOlO ClOll 
\CllOO CllOl ClllO Cllll/ 

Using this we write 

(g 

-eg \ 

Using the special form of the matrix C we have the property (gC^g = hence a staightforward 
calculation shows that 

[eg 1 + 

where 1 = / (g) / and 



(54) 



(55) 



-C^^7C^ 



^(0) . ^(0) ^(0) . ^(1) 
^(0) . ^(1) (^(1) . ^(1) 



(56) 



Here the 4-component vectors C''°^ and C''^^ are just the first and third columns of the matrix 
C of Eq.( |53l) . and the ■ product is defined by Eq. (11851) of the Appendix. 
Due to the special structure of C we also have the property 



(57) 
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resulting in our final form for the coset representative in the Iwasawa gauge 



V 



M3 ® M2 







Ml (g) Mc 











) 




) 



(58) 



We close this section with some important comments. From the particular form of our 
coset representative in the Iwasawa gauge, also reflected in our choice of the matrix ( of 
Eq. (l53|l . we see that the role of the first qubit labelled by io is special. The corresponding 
SL{2, M) action refers to the Ehlers-group. However, our choice of ( also gives special status 
to the second qubit labelled by ii. This is also reflected in the structure of the matrix A of 
Eq.f lSB]) . The 8 components of ( can be regarded as the ones arising from an embedding of 
a three-qubit state sitting inside a four-qubit one having merely 8 nonvanishing amplitudes. 
The grouping of these amlitudes of this three-qubit state into two four-vectors C^^^ and C^^^ 
is based on the special role we have also attached to the second qubit. However, we would 
have chosen any of the remaining two qubits to play this role. This would have resulted 
in another 4 plus 4 split for the 8 nonzero components of (. This freedom for different 
arrangements is related to the triality of so(4, 4) connected to the permutation symmetry 
inherent in the embedded three-qubit system. For more details on this point we refer the 
reader to the Appendix. 

IV. THE LINE ELEMENT ON M3 AS A FOUR-QUBIT INVARIANT. 

The line element on A^3 is given by the formular^ 



ds^ = Tr{ry 



(59) 



where 




(60) 



and the involution compatible with our conventions is 




(61) 
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Using the explicit form for V as given by Eg. (1581) a straightforward calculation gives the 
result for V 

^ ^ 1 /Sg ® /2 + /s ® S2 -g-^ - ^g ^^^^ 
2 \^ c/^^ + ^^g Si ® Jo + Ji ® Sc 

where 

^, = ^(-'"-'^'1, , = 1.2.3 (63) 

Vj \ —dxj dyj 

and 

1 / —dvo —dxn + w\ 1 r ~ ~ T 

So = - , w = -{C'dCi-CidC'). (64) 

X-dxQ + wdyol ^ 

The important part we have not discussed yet is the 4x4 matrix 

^ = (M3 ® M2)dC{Mi ® Mof, (65) 



which by virtue of Eqs.f l471l49|) can be written as a differential form on the symplectic torus 
determined by the Wilson lines based on a four-quhit state 

1^) = (MgOMsOMi (g)Mo)|rfC). (66) 

Recalling our conventions of Eqs. (l50|) . ( 152|) . (!53|) we expect that is depending on the warp 
factor, the NUT potential, the moduli, and the Wilson lines d(^ and d(i. The four-qubit 
state \d() depending only on the Wilson lines clearly determines the entanglement type, 
since is lying on the S'L(2,M)^^ orbit of this state. However, due to the special role of 
our first qubit |\E') is of special kind. Like in Eq.f l53p its nonzero amplitudes when displayed 
in a 4 X 4 array are located in the first and the third columns. An important consequence 
of this is that the NUT potential is not appearing in the explicit form of |\1'). 

We can get a four-qubit state |$) of a more general type after reinterpreting the term 
g'^g + \1/ = [g"^ + '^g)g found in the upper right block of Eq. fl62l) as a superposition 

1$) = (e®£®£®e)|^) + 1^). (67) 

The explicit form of this state is 

1$) = (M3 O M2 Ml ® Mo)\dC) + {Ms (g) M2 O Ml ® Mo)^^ VC) (68) 
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where \() = (e (g) e ® £ (g) £)\0- We see that \() is transforming with respect to the contra- 
gredient action. Using Eq. (l53|) the exphcit form of the transformation \() h-)- \() is 



/-Co Ci o\ 
C2 

Ca 



/o -c° C\ 

-Cs 
-C2 

\Q -Ci -Co/ 



(69) 



i.e. fields with a tilde are transformed into the corresponding ones without a tilde up to 
some crucial signs {(j ^ and ^ ~Ci), and their locations are shifted by one column. 
Using these results the final form of V is 



V 



I2 + S2 -^g 

Si ® Jo + Ji ® So , 



Using Eqs. fl59tl62|) we obtain for the line element the following form 



(70) 



where 



1^1 



dx'j + dy^ {dxo - wf + dy, 



yl 



1^1 



(71) 



e(^I^), I^) = (M3®M2®Mi®Mo)MC)- (72) 

Alternatively we can consider Eq.f TfOl) featuring $ = g'^g + which is the 4x4 version 
of the state |$). Then by virtue of the special structure of the matrix \E' (which is similar 
to the one of Eq. (15H]l ) satisfying '^g'^^g = one gets 



{^\^) = Tr(M/^M/) = ^TT{<^g<^^g) = l£*3.^,^«2.i.^n.',^.o.;,$^^^^^^^^$., ., . (73) 



We see that the term occurring in the expression of the line element has the im- 

mediate interpretation as the norm of a four-qubit state. However, again due to the special 
structure of \dQ which determines the orbit type of |\E') it is natural to give a three-qubit 
reinterpretation as follows. Define 



= (M3 (g)M2 (8)Mi ® /)|rfC), 1^) = (J® /®/® Mo)|^/'). 



(74) 



Then we have 



1^1 



1 

yo 



2 ^ ^-2U 



(75) 
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where by virtue of 

\dC) = dCi3i2ho\hi2iiO) (76) 

WipW"^ can be regarded as the norm squared of a three-qubit state. Let us now recall Eqs. ffT^ - 
(IT^ and One can check that 

i.e. we get back to the usual notation used in the supergravity literature. Notice that unlike 
its usual form the new version as a norm squared is not explicitly SL(2,M.)^^ C S'p(8,M) 
invariant. However, by virtue of Eq.f l73|) we have another interpretation for this term, 
which clearly displays its SL{2,M.)^^ invariance. (The expression is actually the canonical 
quadratic SL{2,'M.)^^ four-qubit invariant, however, the Ehlers S'L(2,M) transformations of 
the form / (g) / / 5* are not preserving the special form of |$).) 
Proceeding further let us define the 8x8 unitary matrix 

(u®u \ I (i l\ 

U= , U^HP=-=\ , (78) 

y u^uj -ij 

where H and P are the Hadamard (discrete Fourier transform) and phase gates known from 
Quantum Information Theory 

We notice that for a = 0, 1, 2, 3 
, i I dza \ 

WEJJ = — dzj = dxj — idyj, dzo = {dxo — w) — idyo, (80) 

y-dz^ J 

with j = 1,2,3. After introducing the right invariant one-forms Cq, = ^dza on the cosets 



[SL{2,R)/SO{2)]^ we can define 



e„ = UT.JA^ = I ° I . (81) 

ea 



Using the unitary matrix of Eq.f l78p we can transform V of Eq. fl62|) to the form 

, 1 / 63 ® /2 + /s ® 62 +"^9 , 

V = VVV^ = - I ^ ^ I ^ ^ I . (82) 

-g^^ - ei (g) Jo + Ji ® Co , 
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Here we have introduced 4f — {U ®U)'^{U answering the new 4-qubit state 



1^) = {H ® H ® H ® H){P ® P ® P ® P){M^® M2® Ml® Mo)\dC,). 



(83) 



Notice that this new 4-qubit state is now on the (71/(2, C)^'* orbit of the one \dC,) due to 
the presence of the matrices U G U{2). Moreover, |^) can also be regarded as the discrete 
Fourier transform of the one {PM^ (g) PM2 PMi ® PMQ)\dC,) incorporating the important 
phase factors e*^ via the phase gates. 

In order to gain some insight into the structure of V we define the 4x4 matrix 



corresponding to the complex four-qubit state 



= (£ (8) £ (8) £ (g) £)|*) + I*). 



(84) 



(85) 



Notice that though this state is now complex it is again of special form since it satisfies the 
reality condition 

1$) = (£(8)£(8)£(8)£)|$). (86) 

In order to understand the structure of |l>) we write its component state |^) in a three- 
qubit-like notation 



1^) = {I®I®I®UMo)\ij), 



(87) 



where again ■^igiano 7^ but ipisi^ni — thanks to the structure similar to that of Cisiiwo- 



Introducing the shorthand 
in 4 X 4 notation we get 



(89) 



'^^000 ^000 ^001 ^001 \ 
i ^010 ^010 ^011 ^011 
v2z/o ^100 ^100 ^101 ^101 

Vv'iio V'lio V'lii V'lii/ 

i.e. the first and the last two columns are the same. Now using the special structure of the 
matrix U ®U ®U one can verify that the following reality conditions hold 



111 



000 > ■0001 



IIOj 



'010 



101 > 



100 



oil- 



(90) 
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As a result of these considerations the matrix $ takes the following form 

(Iq ^\ 



^ ^3 
^ ^2 
V^i £oJ 



(91) 



with 



So = 2^11 



10 



— V'lii, 



1100 



2 7 

—Vno, 

yo 



So = 2^ 



1010 



101) 



S-^ = 2^ 



Olio 



on- 



(92) 
(93) 



After using this result in the expression for V of Eq. fl82|) we arrive at the explicit form 



V 





/ 





62 


ea 





£1 






















63 





-^3 


-^2 









es 








62 





-^2 


-^3 





1 







es 


62 





^0 








^1 


2 




-^1 








-^0 





eo 


ei 














£2 





eo 








ei 











S3 





ei 








eo 




V 


-So 













ei 


eo 


0/ 



(94) 



The line element in terms of these complex quantities is the familiar one of Eq.(l7T 



yl 



1 

yo 



1^1 



2/0 



(95) 



a=0 0=0 ""^ •'^ a=0 

Here = (^|^) is the usual scalar product on C^^ with complex conjugation in the first 

factor. 

It is important to realize that our quantities £a can be written in the familiar form 



£, = V2e'^-''X\MijdC' - dCi), £, = 2zv^|/,e-^// (AfjjdC' - dO) 
in terms of the quantities known from special Kahler geometry. Here 

/l\ 

1 



(96) 



// = e^DiX^ = e^{di + {diK))X^ = — 



Zi - Zi 



e.t.c. 



(97) 
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with = (1, zi, Z2, z^Y , K = — log(?/i?/2?/3) and Mij is defined by Eqs. fll2H13|l . By virtue of 
Eqs.f l92ll93|) tliese quantities are nicely compressed into tlie four-qubit state of Eq. fl83|) . 

Tlie special structure of V of Eq. f l94p reveals yet another way for obtaining a four-qubit 
state. Indeed, V contains precisely 16 nonzero quantities which can be organized to form 
the amplitudes of this new state. In order to motivate our construction of this new state let 
us consider the space of 8 x 8 matrices of the following form 

/O 6 a 0\ 



B 



n 



/-7 - 5 

-a + P 
a- (5 




C 



\ 






V 







A 



7 + 5/ 



c 

c 

\0 a 6 0/ 

I -a -13 

-7 + 5 



\ 



(98) 






7 — 5 




a 



\ 




+ /3 / 



(99) 



where a, /3, 7, 5 G M, and a, b,c,d ^ C. Clearly this 12 real parameter family is complemen- 
tary to the 16 real parameter one characterizing V. In both cases the off diagonal blocks are 
related as C 1— t- —gC^g. In the case of V the off-diagonal block satisfies the reality condition 
$ = g{^)g, and for C this condition is C = —gCg. Both of the matrices V and TZ are 
satisfying Eq.(|38ll hence they are elements of the Lie algebra of 5*0(4,4, C) ~ 5*0(8, C). 

Let us now label the rows and columns of these 8x8 matrices as 0, 1, 2, 3, 4, 5, 6, 7 or in 
binary notation 000, 001, 010, Oil, 100, 1010, 110, 111. Now we employ the following permu- 
tation to the rows and columns 



(0, 1, 2, 3, 4, 5, 6, 7) ^ (7, 1, 2, 4, 3, 5, 6, 0), (100) 

(000, 001, 010, Oil, 100, 101, 110, 111) ^ (111, 001, 010, 100, Oil, 101, 110, 000). (101) 

The binary notation is instructive since it clearly shows that after applying the permutation 
we get two 4 element blocks labelled by numbers containing an even number of zeros for the 
first block and an odd number of zeros in the second. (Another mnemonic: the numbers 
1, 2, 4 are the quadratic residues modulo 7 and the ones 3, 5, 6 are the quadratic nonresidues.) 
Now it is easy to check that our fundamental matrix G of Eq. ( 139|) is invariant under this 
permutation. 
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Applying this permutation to the matrix 1Z yields the one 





/ 

'7 + ^ 


c 


d 






















c 


-7 + (^ 





d 



















d 


7 — (5 c 

















2 


U 







c 



—7 — 5 





a + P 


U U 

a h 


U 





















a 




-a + (3 


b 




















h 




a-P 


a 
























b a — 


a — 






This matrix contains the two 4x4 blocks in its block 


diagonal part 












u> 
















5 c<t), 


d = 


d2 


, c = 






d — di + id2, 


c 


= Ci + iC2 








VJ 




vv 
















fh] 














§ ao"), 


b = 


h 








b = bi + ib2, 


a 


= ai + ia2 








VI 















(102) 



The same permutation acting on V results in the new form 




(105) 



^Annnn An 



A 



An 



Anni 1 ^ 



^0000 1^0001 ■'1^0010 ■'1^0011 
^0100 Aqioi Aqiio Aqiii 

Alooo Alool Alolo Aloii 
\Aiioo Aiioi Aiiio Aiiiiy 
Now we define a new four-qubit state 



( -Sq -Co -ei 

62 S2 Si 63 

63 82. 82 62 

\-8i -ei -Co -80 J 



(106) 



|A) 



A. 



03020100). 



(107) 



03,02,01,00=0,1 



Looking at the structure of our matrix TZ' it is clear that it defines an infinitesimal SU{2y 
action on our state based on the 4x4 complex matrix A related to the decomposition 
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so(8,C) = [s/(2,C)]^© (2,2,2,2), 

and the embedding of su{2) in sl{2, C). 

It is important to realize that after the transformation 

fo -b -a \ 

d -c 



(108) 



c -d 
\0 a b / 



(109) 



with C having the property C = gC'g and also the one of Eq. fllOip the matrix replacing 
Eq.( ll02p will contain the diagonal blocks 



1 



^dr (g) / + / (g) cr). 



1 



br (g) / + / (g) ar). 



Here the matrices 



r2 ^ 



rs = -as, 



(110) 



;iii) 



satisfy the commutation relations of an sm(1, 1) subalgebra of sl{2, C). 

Recall that one of our qubits (i.e. the first one labelled by ao) is still special. According 
to Eq.f llOGp transformations of the form I ® I ® I ® S, S = e^^" G SL{2, C) acting on this 
qubit relate the first column with the second, and the third with the fourth. Hence these 
transformations relate Sq to Cq and the Sj to the ej with the same index j. 

Notice that the SU{2) subgroup of this S'L(2, C) is just the R-symmetry group arising 
from the restricted holonomy group, of the quaternionic-Kahler space which is the analyti- 
cally continued version of our para-quaternionic M.^- This holonomy implies^Is^i^ that the 
complexified tangent bundle of that space splits locally as W (g> V where W and V are vector 
bundles of dimension 8 and 2. In our case a similar split exists where the former space corre- 
sponds to the three-qubit part of our A labelled by the indices 03, 02, ai and the latter to the 
special qubit labelled by Oq. Indeed our transformations of Eqs. (!78|) and (110 ip correspond 
to a change of basis in T^Ai^ similar to the usual one rendering the " quaternionic vierbein" 
covariantly constant with respect to the spin connection^!. 

Eqs.f ll05lll07p are of central importance for our considerations of the following sections. 
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They define a complex four-qubit state satisfying tlie reality condition 

|A) = ((Ti ®(Ti ® di ® (Ti)|A), (Ti = , (112) 

V V 

wliere cxi is tlie bit flip gate of Quantum Information Theory. It is straightforward to check 
that the subgroup of transformations of the group SL{2,C)^^ leaving invariant this reality 
condition is SU{1, 1)^^ i.e. precisely those transformations as described by Eqs.f ll09|) - flllll) . 
Hence in the notation used in quantum information the admissible transformations are of 
the form 

\A) ^ iS^S2(S)Si^So)\A), Ss,S2,Si,Soe SUil,l). (113) 

In what follows our basic concern will be a study of quantities invariant under the larger 
group of transformations i.e. S'L(2,C)^^. Such invariants are clearly also SU{1, 1)^^ ones. 
It is known that the number of such algebraically independent invariants is four— We 
have a quadratic, two quartic, and one sextic invariant. The structure and geometry of such 
invariants has been investigated^!^. Here in closing this section we just observe that the 
quadratic four-qubit invariant^ for our state |A) is precisely the line element dsj^^ i.e. 

1 3 

r/s^ = ^aga' aaa' aia' aoaj, A \ , , , , = \^ (p e -E E\ fll4) 

This formula flrst appeared in the paper of Bossard et.al.— Here we have also clarifled its 
intimate connection to four-qubit systems. We also remark that the quadratic invariant is 
also a permutation invariant. However, from the physical point of view the special role we 
have attached to the first qubit obviously breaks this permutation invariance. 

V. CONSERVED CHARGES 

The 3D duality group acts isometrically on our A^3 by right multiplication and yields a 
conserved Noether charge^>2^>^ 

Q = V-'VV = I I . (115) 



The explicit expression of Q is given by 



-eg l + ^Aj \ gdC^ + dCg p,®I, + h®p,) y^g l + ^A, 
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where 



dz, 



= Re I — I 1(2:^1)1, Za = Xa-iya, a = 0,1,2,3 (117) 



and 



rfC = (ATg® Ar2)dC(A^i® A^o), Ar„ = MjM,, a = 0,1, 2, 3 (118) 

where for dzj,j = 1, 2, 3 and dzo we used the definitions of Eq.flHOl). 

Now we are interested in the conserved electric and magnetic charges coming from the 
first and third column of Qi2- (This part has the same structure as the matrix ( of Eq.( l53|) ). 
Since the matrix (g has vanishing first and third column, any 4x4 matrix multiplied by (g 
from the right also has this property. Hence terms having this structure will not contribute 
to the relevant part of Qi2- Using A = —('^gCg the relevant part of Q12 is 

2[gi2]relevant = "(iVs ® N2)dCiN^ ® /) + \ {c ® £)((£ /)• (119) 

Uo Ho 
On the other hand let us look at the conserved quantity 

k ^ ^Tr(J, ® E,)Q22 = i^^^^^. (120) 



Since the line element is related to the Lagrangian and the Lagrangian according to Eq. flTT]) 
contains a term ((io — ui)'^ + l/o^)/4?/o hence = Pa = §^ = k i.e. our quantity is just 
the NUT charge. These considerations show that the relevant part T of Q12 in four-qubit 
notation is 

|r) = ^iNs^N2®N,^I)\dC)-Pxoie^e(g)e(g)I)\C). (121) 

Eq. fll2ip comprises 8 conserved quantities represented as the 8 nonzero amplitudes of a 
four-qubit state. Note, that our formula also contains the NUT charge. In a three-qubit-like 
notation we can alternatively write this as 

ID = ^e-2^(Ar ® I)\dO - p.(e ® /)|C), (122) 

with 

Af = Ns^N2^Ni, e = e(g)e(g)e. (123) 

Let us see how this set of conserved quantities is related to the momenta p^i = ^ and 
PCi ~ f^- ^ calculation based on the Lagrangian related to the line element Eq. fl71l) shows 
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that these quantities can be also organized into a state 

be) = -^e-^^lA/" ® /)MC) + y (e ® /)|C). 



Hence 



|r) = -|(e®/)|C)-bc)- 



Let us now introduce the new quantity 



Pa 



(124) 



(125) 



(126) 



(127) 



in) = bc)-Y(e®^)IC). 

Notice that after writing out the 8 amphtudes exphcitly we get 

7-) Per T „ Per 

^c^ = PC - yC/, = Po + yC , 

in accordance with Eq. (4.14) of Bossard et.aL— (The a used by them is different by a factor 
of 2). Now one can verify that the Hamiltonian is 

3 

H = Y.y'(Pl +pI) - yoiPcl^-' ® I\Pc)- (128) 



0=0 



From Eqs.f lT25lfT26|) we have 



|f) = |r)+p.(e®/)|C) = -|Pc) 



(129) 



arriving at an alternative expression as 



H = J2yl(pl+pl)-^"''im-'^m 



a=0 



in a three-qubit-like notation. 

Let us now parametrize F in terms of the electric and magnetic charges as 

/ / -p^ o\ 



1 

7! 



gg 
-p^ g2 
\qi go 0/ 

where the rows and columns of this matrix are related to the four-qubit labels as 
defining the four-qubit state |F). This parametrization in the conventional languag 
to 



(130) 



:i311 



1 



V2 



p^ = ^p.(^ + Pcj ' -^1^ = \paO - PC 



in Eq.(j53!) 

;e amounts 

(132) 
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Notice that the origin of the factors of \/2s appearing in Eqs. fll3mi32]) can be traced back to 
the fact^i^ that the electric and magnetic charges should be proportional to the generators 
\/2Epi and y/2Eqj. For vanishing NUT charge Pa- = we get 

c^^Vbh = e2^(r|A/'-i ® /|r) (133) 

which can be checked to yield the usual expression for the Vbh black hole potential. 

Let us now rewrite our expression of Eq.f l66|) for in terms of our conserved quantities. 
First by using Eqs. (11221) and (I129p we express \d() in terms of the charges as 

|dC) = 2e2^(Ar-i® J)|f), (134) 

to arrive at the expression 

1^) = 2e^{Mj ® Mj ® ® /)"Nf ). (135) 
We can further transform this to obtain |^) = (W (g) W (g) W (g) W)|^) of Eq.(IH3]) 

1^) = 2e^^{I / ® / ® WMo)(V ® V ® V 8) /)(53 ® ^2 ® 5i ® /)|7), (136) 
where we have used that 

UM^'^ = VSas. (137) 

Here the matrices V and Sj, j = 1, 2, 3 are the ones of Eq.(l2]) discussed in the Introduction 
and 

I7) = (a3®(T3®(T3®/)|f). (138) 

Clearly for A; = i.e. vanishing NUT charge I7) is just the phase-flipped version of |r) of 
Eq.( ]13ip . This I7) reinterpreted as a three-qubit state is just the charge state mentioned in 
Eq.dl]). Now notice that we have 



1 * 



UMo = -j=e-^ ^ ' ' ^^^^^ 

where the terms in the second column are not needed since for the four-qubit state |r) we 
have as usual t'j-jj-jjii = 0, a structure that dates back to the similar one of \() and |r). 
Now obviously ^l/jgjaiii = "^jsjijiO and either of them can be reinterpreted as the ones of a 
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three-qubit state. (See also Eq. (189|) in this respect.) For this three-qubit projection we 
have 

1^)3 = iV2\x) = iV2e^{V ® V ® V)(^3 ® ^2 ® ^i)|7), (140) 
where by virtue of Eqs. f ll38p and (11291) 

I7) = l7)+P<x(^i®^i®^i)|C)- (141) 

Eqs. fll40p and (I14ip clearly show how our state \x) of Eq. ([7]) as a special case of \x) 
is embedded in a four-qubit state |^). Moreover, in achieving this we managed to present 
a generalization also valid in the case of nonvanishing NUT charge. The important new 
property to be noted here is that unlike I7) which is constant the one I7) is depending on 
r = i via the Wilson lines and (j. Notice that we also have = 'ie^^VsH = 4||x||3) 

i.e. the Black Hole potential is just the norm of a three-qubit state.->^. 

For later use for static spherically symmetric solutions let us write out eplicitly the 
quantities £a of Eq.( P^ - (l95]) in terms of the three-qubit state \x) of Eqs. ([7]) and ([T]) 

£0 = iVSxni, £\ = iVSxuo, ^2 = iVSxwi, £3 = iVSxou, (142) 



£o = w8xooo, i^i=w8xooi, ^^2=«v8xoio, i^s = w 8x100 • (143) 

For the more general stationary solutions we have to use \x) as given by Eq. (11401) . 
In closing this section let us also calculate the conserved quantity 

m = ^TT{h®Ho)Q22. (144) 

Writing out explicitly Q22 using Eq. (lll6l) we notice that many terms end with the matrix 
(g. Using the cyclic property of the trace these terms in some cases result in ones begining 
with Cg{Ii (8> Ho)(^g which is vanishing. Employing Eq. (ll32p and the definition = ^ 
the result of these considerations will be just two nonvanishing terms yielding the final result 



^ = ^(r|C) + xop^o + yoPyo = -^iC^Pc' + OPq) + (^P<T + U. (145) 
which is the ADM mass of the black hole^"— . 
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VI. BLACK HOLE SOLUTIONS AS ENTANGLED SYSTEMS 



A. BPS solutions 

Let us consider our four-qubit state |A) of Eqs. fllOSp - (11071) . In this section we would like 
to investigate issues of separability for this state. In particular in this subsection we will be 
interested in the sufficient and necessary condition for the separability of the first qubit, i.e. 
the one which is labelled by ao in Eq. (11071) . From our previous considerations it is clear that 
this qubit is the one of special status, i.e. it is the one transforming as a doublet under the 
i?-symmetry. 

In QIT terms separability of this qubit from the rest is equivalent to the condition that 
the (unnormalized) 2x2 reduced density matrix gi = Tri|A)(A| represents a pure stated. 
This density matrix is of the form 

/(Ao|Ao) (Ao|Ai)\ 

^1=1 , , , , I > (AaolAaJ,) = > K3a2aiaoK3a2aia'^- (146) 

\^(Ai|Ao) (Ai|Ai)y a3,a2,ai=0,l 

This is a pure state i.e. a projector if and only if Detf^i = 0. Equivalently this condition is 
satisfied iff the two 8 component vectors Aagaja^o and Aagaaaii are proportional, i.e. AagajaiO = 
AA(i3(i2aii- By virtue of the reality condition of Eq. (lll2p we also have the constraint |A| = 1. 
Using the definitions in Eq. (ll06p this means that 

Sq = Aeo, Sj = Xej, |A| = 1. (147) 

The first consequence of these considerations is that for the state A the quadratic invariant 
Ji of the Appendix (see Eq. (11890 ) which is related to ds"^ of Eq. (lll4p is vanishing. Moreover, 
since the first column of the matrix A of Eq. (ll06p is proportional to the second and the third 
one is proportional to the fourth, the invariant J4 which is according to Eq. (ll9ip just the 
derminant of A is also vanishing. A straightforward calculation based on Eqs. (ll90p and 
(I193p shows that the remaining two algebraically independent invariants I3 and I2 are also 
vanishing. 

Now the 4x4 matrix Q = A^gAg satisfies Eq. (ll97p of the Appendix, so it follows 
that = 0. This implies that the matrix TZ\ of Eq. (I202p is nilpotent. According to 
the terminology of Ref.— states with the property that their associated 8x8 matrix 7?.a 
is nilpotent are called nilpotent states. Hence our separable state is a (trivial) example of 
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a nilpotent 4-qubit state. Moreover, since T^a is just 2V' of Eq.f llOSp and this matrix is 
unitarily related to the matrix Q of Eq. flllSp of conserved charges, it follows that Q is also 
nilpotent. 

In order to link these considerations to the usual static, extremal spherically symmetric 
BPS black hole solutions we choose A as 




Note that for static solutions the NUT charge is zero, hence = and = = dcf), 
i.e. = Co. 

Now in the language of supergravity the above discussed condition on separability is just 
the usual one on the existence of Killing spinors2ii22 expressed in terms of the quaternionic 
vierbein 




(149) 



and Eqs. fll47p . (l96l) and ( l3Ti) . give rise to the attractor flow equationai^i^I 

U = -e^\Z\, = -2e^G^'^d^\Z\. (150) 

As it is well-known these first order equations imply that the corresponding second order 
equations of Eq. fl33p hold. Moreover, by virtue of the vanishing of the invariant Ji the 
constraint of Eq.( l28|) is also satisfied hence the solution is extremal. 

From this analysis we have learnt that the condition of separability for the first qubit for 
the 4-qubit state |A) taken together with the special choice of Eq. fll48p yields the first order 
attractor fiow equations. Moreover, we have seen that in this case |A) is a nilpotent state. 
This property of |A) is related to the well-known nilpotency of the Noether charge2Ii2^i2i Q_ 
Notice however, that our approach does not directly yield the order of nilpotency of Q for 
BPS solutions which is three^. 



B. Non-BPS solutions with vanishing central charge 

Let us discuss the separability properties of |A) associated with the remaining qubits not 
playing any distinguished role. Here we chose to consider separability of the 4th qubit. An 
argument similar to the one as given in the previous subsection shows that the sufficient and 
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necessary condition of separability for this qubit is that the first row is proportional to the 
third and the second is proportional to the fourth. Due to the reality condition we again 
have |A| = 1 and we get 

£q = -Aes, Si = -Xe2, £2 = -Ae2, S3 = -Xe^. (151) 

Using the definitions of Eq. (!96l) these conditions take the explicit form 

^ = Xe^Zs, ^ = Xe^Z2, ^ = Ae^Zi, ^ = -Xe^Z, (152) 

yo yi y2 ys 

where Zj = —2iyjDjZ with Dj as given by Eq. (!32|) . Now for static solutions we again have 
no twist potential i.e. Xq = hence by choosing 

X = -i 

we get 

U = -e^\Zs\, = -2e^G^^d^\Z3\. (154) 

These expressions show that demanding separability for the fourth qubit taken together 
with the choice of Eq. (11531) yields the first order equations characterizing attractors with 
vanishing central charge^. 

Clearly similar considerations apply for issues of separability for the second and third 
qubits. The result will be similar sets of equations with IZ3I replaced by \Zi\ and \Z2\. This 
amounts to taking different forms for the fake superpotential^. 

Note that the value for the four-qubit invariant Ji is related to the extremality parame- 
ter. Unlike the other three algebraically independent invariants, this is also a permutation 
invariant. Of course the value of Ji is zero for both BPS and non-BPS solutions with van- 
ishing central charge, expressing the fact that our solutions are extremal. Moreover, for all 
of our non-BPS solutions some rows or columns of the 4x4 matrix are proportional, hence 
the invariant J4 is zero as well. Calculations show that the remaining invariants I2 and I3 
also give zero, hence our considerations on the nilpotency of |A) familiar from the previous 
subsection still apply. 

In closing this subsection we note that the conditions for separability can be written in 
the familiar formal of Eq. (11491) with the label of e"" is Oq for BPS, aj, j = 1, 2, 3 for non-BPS 
solutions with vanishing central charge. Of course A should be modified accordingly. 
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C. Non-BPS seed solutions 



From the previous subsections it is obvious that the condition of extremahty related to 
the vanishing of the invariant Ji can be satisfied in a number of different ways. Exphcitly 
the relevant equation to be satisfied is 

3 3 

«=0 a=0 

For static solutions we have already remarked that cq = eo, hence for BPS solutions Eqs. 
f ll47p - fll48p can be written in the form = Ae^, i.e. £a is related to via a special element 
of U{A) containing merely phase factors A in its diagonal. In the case of non-BPS solutions 
with vanishing central charge Eqs.f ll51| - fll53p of the previous subsection can be written in 
a similar way in terms of another element of U (4) 
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Similarly the basic equations of the remaining two cases of the previous subsection can 
be expressed in terms of similar unitaries. These unitaries are just permutation matrices 
combined with phase factors and their conjugates. As we have shown this structure is related 
to the separability of one of the qubits of the state |A). In simple terms this means that 
some of the rows or columns of the 4x4 matrix A corresponding to |A) are proportional to 
each other. 

In order to obtain states |A) which are entangled and at the same time give rise to static 
spherically symmetric non-BPS black hole solutions with non-vanishing central charge we 
have to experiment with elements of f/(4) of more general type. 

Let us consider the following choice 
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Due to the unitarity of the relevant matix the condition of extremality is satisfied, moreover 
obviously none of the qubits can be separated from the rest. However, apart from satisfying 

31 



Eq.f llSSp or equivalently Eq. (l28|) we still have to satisfy the equations of motion i.e. Eq.( l33|) . 
In the following we show that the choice of Eq. fll57p indeed gives rise to a solution of the 
latter equations namely the non-BPS seed solution^^. Clearly apart from characterizing the 
seed solution in a nice and compact way Eq. fll57p also serves as a mnemonic for the structure 
of the corresponding entangled 4-qubit state |A) of Eqs. fll05l) - fll07p . 

In order to reveal the structure of the seed solution for special non-BPS charge con- 
figurations we recall Eqs. fll42p -f ll43p and ([7]) and employ a discrete Fourier (Hadamard) 
transformation to \x) as 



\x{T)) = {H^H^H)\x{r)), H = —\ . (158) 




The amplitudes of this state are 

\/2?/i2/22/3Xooo = -ie" ymy-iP^ , (159) 
V^^ymmXiw = ie^?/i(x2X3p° - X2P^ - xsp"^ + qi) (160) 



x/^ymysXwi = ie^y2{xix^p^ - xip^ - xsp^ + ^2), (161) 



V^ymmXou = ie^yz{xiX2P^ - xip^ - X2P^ + gs), (162) 



V^mymXui = e^{XiX2X3P^ - X2X3P^ - XiX^p^ - XiX2P^ + XiQi + X2q2 + 2:3^3 + go), (163) 

V^ymrnXooi = e^2/22/3(p^ - xip^), (164) 

V^yiymXow = e^ymip"^ - X2P°), (165) 

V^yiymXooi = e^yiy2{p^ - X3p°). (166) 

Now one can check that Eq. fll57p can be expressed in terms of these quantities as 

_ixo ^ _ixi ^ _ix2 ~ _ix3 
2^/0 2 2/1 2 2/2 2 2/3 

1/2/0 2/1 2/2 2/3\ ~ 1 / 2/0 , 2/1 2/2 2/3\ 

Xiu = 7 , Xooi = T + (168) 

4 V2/0 2/1 2/2 2/3/ 4: 
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2/1 



1 / 2/0 2/1,2/2 2/3 \ ~ 1 / yu tfi tfz , 1 /1cn^ 

Xolo = 7 + , Xloo = 7 + — • 169 

4 V 2/0 2/1 2/2 2/3/ 4 V 2/0 2/i 2/2 2/3/ 

For static solutions we have vanishing NUT charge i.e. xq = hence the first of these 
equations reads as xooo = which by virtue of Eq. fll59p means that p° = 0. Hence our 
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candidate for a non-BPS solution should have only seven nonvanishing Fourier amplitudes 
and no D6 brane charges (in the type IIA duality frame). 
Let us now introduce the notation 

2/o = e<^°, % = e<^^ /3 = f/-^(0i + 02 + 03), a, = U+^(f),, (170) 

with and j = 1,2,3 (recall also that according to Eg. (1501) now = 0o = 2f/.) Now our 
equations take the form 

X ' z z % 

Xiii = Xiio = ^e~'>"xi, xioi = ^e~'>'^i2, Xou = ^e~'^''x3, (171) 

Xooi = ^(«i - «2 - as), Xoio = ^("2 - "3 - ai), Xioo = ^(«3 - «i - «2)- (172) 

Now using Eqs.f ll60| - fll66p with the further charge constraints g-,- = , go < 0, and 
p^,p'^,p^ > one can see that the equations are precisely the ones found in the Appendix of 
the paper of Gimon et.al.^^ characterizing the seed solutions for the DO — D4 system. 

We remark in closing that one can verify by an explicit calculation that all of the four 
algebraically independent four-qubit invariants Ik,k = 1,2,3,4 are vanishing. This means 
that the corresponding matrix Q of conserved charges is nilpotent. Hence in the teminology 
of four-qubit entanglement we obtained the result: the relevant state |A), is a nilpotent one. 
However, unlike in the previous cases now neither of the qubits can be separated from the 
rest, hence |A) is also an entangled state. 

Notice however, that neither the order of nilpotency nor the particular entanglement type 
follows from our simple considerations. It would be interesting to extend our analysis and 
identify the particular entanglement class to which |A) belongs case by case. It is important 
to realize in this respect, that in our simplified considerations we have merely used complex 
four qubit states and the corresponding S'L(2,C)^^ invariants. However, we must recall 
that our state |A) also have to satisfy the reality condition of Eq. (11121) . The result of 
the implementation of this constraint is that in the black hole context we have to classify 
orbits under the group SU{1, 1)^^ which is merely a subgroup of the full group of admissible 
local operations^. Hence a full entanglement based understanding of black hole solutions 
in the STU model should rely on the classification of entanglement types of real four qubit 
states defined by Eq. (I112p . This classification should be founded on a study of SU{1, 1)^*^ 
invariants. Clearly this reformulation would relate the known classification^! of satic black 
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hole solutions in the STU model in terms of nilpotent orbits to a similar one based on 
entanglement classes of the relevant real four qubit states. 

VII. CONCLUSIONS 

In this paper we managed to understand the structure of extremal stationary spherically 
symmetric black hole solutions in the STU model of D = 4, = 2 supergravity in terms of 
four-qubit systems. Our analysis extended the results obtained in our previous papers based 
on three qubit systems^i^ii^. The basic idea facilitating this 4-qubit description was the fact 
that stationary solutions in D = A supergravity can be elegantly described by dimensional 
reduction along the time direction^. In this picture stationary solutions can be identified 
as solutions to a D = 3 non-linear sigma model with target space being a symmetric space 
G/H with H non-compact. The group G extends the global symmetry group G4 of D = 4 
supergravity, by also incorporating the Ehlers SL{2, M). In our specific case the N = 2 STU 
model can be regarded as a consistent truncation of maximal N = 8, D = 4 supergravity 
with G4 = Ejij), truncating to G^ = SL{2,'R)^^. Timelike reduction then yields the coset 
Es(8)/ SO* (16), or in the case of the STU truncation the one M3 = 50(4, 4)/5L(2, M)><^. 
We have shown that the four copies of SL(2, M)s occurring in this coset can be reinterpreted 
as the group of local operations acting on four qubits subject to special reality constraints. 
Here the fourth qubit which accounts for the Ehlers group played a special role. 

The central object of our considerations was the complex A-quhit state |A) of Eqs. (1106P 
and (11071) . also satisfying the reality condition Eq. (lll2p . The amplitudes of this state of 
odd parity contain the right invariant one-forms Cq, a = 0, 1,2,3 defined by Eqs. (l80l) - (!8T!) . 
On the other hand the 8 amplitudes of even parity are just the 8 amplitudes of the 3-qubit 
state well-known from previous studies concerning the black hole qubit correspondence. 
According to Eqs. (19 3 p and ( PB]) these amplitudes are related to well-known quantities of 
special geometry. We have shown that the state |A) is connected to the line element on 
Ais via Eqs. ( [59|) . (llOSp and (HMD . We also realized that this expression for the line 
element is minus the quadratic 4-qubit SL{2,C) invariant Ji of Eq. (ll89p . After expressing 
the 8 amplitudes of the embedded 3-qubit state in terms of the conserved electric, magnetic 
and NUT charges as in Eqs. (ll40p . (I14ip this invariant also has the physical interpretation 
as the BPS parameter—. (For nonrotating solutions this parameter is just the extremality 
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parameter.) 

We clarified the relationship between the warp factor, moduli and charge dependent 3- 
qubit state of Eq.([T]), and ([7]) used in previous studies^'^'i^ and the 4-qubit one |A). Our 
considerations enabled a formal generalization for this state (see Eq. fll4ip ) also valid for 
nonvanishing NUT charge. Notice that for general stationary solutions the entanglement 
type of this state (i.e. the value of the three-tangle^ii^iiS) is also depending on the Wilson 
lines and (j. This is in sharp contrast to the static case where the entanglement type is 
merely depending on the conserved electric and magnetic charges. 

Note that one of the qubits of the state |A) was special. We have seen that the special 
status of this qubit is related to the i?-symmetry group arising from the resticted holon- 
omy group of the para quaternionic Kahler space Ais- We realized that our special set of 
transformations, based on Hadamard and phase gates and permutations, resulting in the 
explicit form for |A) correspond to the basis transformations similar to the ones rendering 
the quaternionic vierbein covariantly constant with respect to the spin connection^!. 

The separability properties of this special qubit are related to the solution being BPS 
or non-BPS. We demonstrated within our formalism the observation of Bergshoeff et.al.— 
that static, extremal BPS and non-BPS-solutions with vanishing central charge^ correspond 
to states for which one of the qubits is separable from the rest. On the other hand using the 
non-BPS seed solutionis for nonvanishing central charge we have shown that |A) in this case 
is entangled. We revealed a connection between the classification of nilpotent states within 
the realm of quantum information theory and the similar classification of nilpotent orbits. 
The details of this connection should be explored further. 

It is amusing to see that nonextremal solutions should correspond to states which are 
semisimple^^ . Since nilpotent states are rather exceptional among the 4-qubit ones, semisim- 
ple states are the ones that represent genuine 4-qubit entanglement. According to our Ap- 
pendix for such states at least one of the algebraically independent invariants (namely Ji 
related to the extremality parameter) are non-vanishing. Such states with special entangle- 
ment properties should correspond to nonextremal solutions of special kind. 

Notice also in this respect that in the paper of Chemissany et.al.— dealing with the full 
integration of black hole solutions in symmetric supergravity solutions the authors notice 
that in their Lax pair approach exactly non-extremal solutions are easier to describe than 
extremal ones. Such solutions correspond to diagonalizahle initial conditions in terms of the 
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Lax matrix. On the other hand they note that the initial conditions summarized in nondi- 
agonalizahle Lax matrices correpond to extremal BPS and non-BPS solutions, however such 
Lax matrices represent a subset of measure zero within the space of Lax matrices. Within 
the context of the STU-model clearly semisimple states should correspond to diagonaliz- 
able Lax matrices, and nilpotent states to nondiagonalizable ones. Hence there should be 
a correspondence between giving the Lax matrix at some initial time and specifying the 
entanglement properties of the corresponding entangled 4-qubit state. 

Notice in particular the highly symmetrical nature of the genuine entangled 4-qubit class 
(see Eq. (l2U4p of the Appendix and the structure of its invariants.) This state is the 4-qubit 
analogue of the famous GHZ state^S familiar from 3-qubit entanglement. Notice that for 
choosing a,b,c,d G M this state automatically satisfies the reality condition of Eq. flll2p . 
Acting on this state with SU{1, 1)^^ transformations preserving this reality condition and 
also the values af the algebraically independent invariants results in a state |A) containing 
16 real parameters. Using this parametrization and the black hole qubit correspondence it 
would be amusing to find a corresponding highly symmetrical non-extremal solution. 

Recall also the classification of black hole solutions in the STU model in terms of three 
qubit entanglement classes as given by Kallosh and Linde^. In this paper the authors noticed 
a similarity between the classification of complex three qubit states'*^ and the corresponding 
classification of small and large black holes in the STU model related to real three qubit ones. 
In the light of our results we might substantially generalize this interesting result. Indeed, 
by embedding the usual three-qubit picture into the four qubit one as described here, we also 
have the possibility to include such notions as BPS and non-BPS , extremal and non extremal 
solutions into an entanglement based picture. As we have seen the extremality parameter is 
related to the quadratic four qubit invariant. Extremal black holes are characterized by the 
vanishing of this quantity. Though the remaining four qubit invariants are all vanishing for 
the known extremal solutions, but such solutions are still distinguished by their entanglement 
properties. For BPS and non-BPS solutions with vanishing central charge one of the qubits 
is separable from the rest, and for the Z ^ case none of the qubits is separable. Since 
the states describing such solutions are real (i.e. they are satisfying the reality condition of 
Eq.f lll2p ) in order to classify their orbit structure we also have to include some additional 
SU{1, 1)^^ invariants. We know that the classes in question for extremal solutions are just 
the nilpoten orbits classified in the paper of Bossard et. al.— , hence these classes might 
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be distinguished by additional SU{1, 1)^^ invariants whose physical meaning is still to be 
clarified. On the other hand the three qubit part of our four qubit states is classified by 
the value of Cayleys hyperdeterminant i.e. the three-tangle^^. For small black holes this 
invariant is vanishing and for large ones it is nonzero and its value is proportional to the 
black hole entropy. 

Finally notice that we have deliberately emphasized the possibility for reformulating the 
well-known results of the STU model in a suggestive permutation invariant language (see 
e.g. Eqs. (!7T]- (!72|) ). Though instructive, this language is deceptive due to the special role 
we have attached to our first qubit via the use of the Ehlers SL{2,M.). However, since the 
quadratic invariant of Eq. flll4p related to the line element and the extremality parameter is 
a permutation invariant quantity one might speculate whether there is a further possibility 
for embedding the STU model into an even greater picture where permutation symmetry is 
manifest. In this respect an exciting possibility is to find the physical relevance (if any) of 
the permutation invariant quantity of Eq. (l200p i.e. the four qubit generalization of Cayley's 
hyperdeterminant (the "four-tangle"). 

VIII. APPENDIX 

A four qubit state can be written in the form 

\-^) = ^ ^i3i2nioK3«2«i«o), \i3i2i1io) = \H)^\i2)<S)\ii)(S)\io) eVi^V2^Vi(g)Vo, (173) 
4312*1*0=0,1 

where V32,i,o = C^- Let the subgroup of stochastic local operations and classical 
communication^ representing admissible fourpartite protocols be S'L(2,C)^^ acting on 
as 

1^) ^ (53®52®5i®5o)|^), S^eSL{2,C), a = 0,1, 2, 3. (174) 

Our aim in this appendix is to give a unified description of four-qubit states taken together 
with their SLOCC transformations and their associated invariants. As we will see states 
and transformations taken together can be described in a unified manner using the group 
50(4,4, C). 
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Let us introduce the 2x2 matrices 
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Ell 




(175) 



Then we arrange the 16 complex amphtudes appearing in '^i^i^i^i^ in a 4 x 4 matrix in three 
different ways 
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(176) 



(177) 



L>3(*) 



*3,*2,n,io=0,l 



«3«2«lio-^«2i3 ® Ei^ig 



(178) 



X W 
Y Z 

where the 2x2 matrices X, Y, W, Z are introduced merely to illustrate the block structure of 
the relevant matrices. The first matrix is obtained by arranging the components of X, Y, W, Z 
as the first, second, third and fourth rows. Notice that the arrangement Di of Eg. (11761) is 
our one of Eg. (1411) . Clearly changing Di to D2 or to D3 corresponds to the two generators 
of the permutation group ^3 acting on the last three gubits. Applying such permutations to 
the gubits corresponds to a similar permutation of the entries S3, S2, Si of Eg. (l40|) resulting 
in the matrices V{s3, Si, S2, Sq] D2), and V{s2, S3, Si, sq; D^). These permutations give rise 
to alternative forms for the matrix exponentials of Eg. (15511561) with special roles attached to 
the second and the third gubit respectively. 

Our matrix V = V{s3, S2, Si, sq] Di) of Eg. lHOll in the parametrization used in Bossard 
et.al.— takes the following form 
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It can be checked that the matrix SVS'^ where 
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is just the one used in Eq.(4.6) of that paper. This matrix also relates our matrix G of 
Eq. fl39l) to the usual 5*0(4,4) invariant one 




(181) 



Relating the upper right block of Eg. (11791) to the 4x4 matrix Dg of Eq. fHOl) shows that in 
this parametrization 
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Epi Eq^ FpO ) 



which justifies our parametrization of Cj^Eq^ + C,iE„i used in Eq.([5 
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Let us discuss now the structure of four-qubit 5'L(2,C)^^ invariants^"^'^'^. The num- 
ber of algebraically independent four-qubit invariants is four. We have one quadratic, two 
quartic, and one sextic invariant. In our recent paper— we investigated the structure of 
these invariants in the special frame where two of our qubits played a distinguished role. 
Clearly this is the case in the black hole context, since one of the special qubits is associated 
with the Ehlers-group and the choice of the other is just a matter of convention related to 
the special choice Di, D2 or of Eqs. (11 7611178 1) . 

To an arbitrary state 

I A) = ^ K3i2iiio\^3i2iiio) , (183) 



«3*2n«0=0,l 



we can associate the 4x4 matrix 
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\Aiioo Aiioi Aiiio Aiiiiy 

or four four- vectors. The splitting of the amplitudes of |A) into four four- vectors reflects 
our special choice for the distinguished qubits compatible with our conventions. Now we 
introduce on the vector space ~ x corresponding to the third and fourth qubit 
a symmetric bilinear form g : x — )■ C with matrix representation: g = e ® e. This 
means that we have an S'L(2,C)^^ invariant quantity with the explicit form 



g{A, B) = gafsA'^B^ = A^B" = A-B = A^B^ - A^B^ - A^B'^ + A^BK 
We can also introduce a dual four-qubit state 

\^) = ^i3i2iiio\^3'i2'il'io) 
«3«2nio=0,l 

with the associated matrix 

/ '^oooo -^0001 -^0010 ^ooii\ / 
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\Aiioo Aiioi Aiiio Ami/ 



(185) 
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where 



a° = e^^^^BpC^Ds, b^ = e'^^^^A^C^Ds = e'^^^^A^BpDs d^ = e^^^^A^B^C^. (188) 
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Here e^^^^ = +1, and indices are lowered by the matrix of g. Notice that the amphtudes 
of the dual four-qubit state are cubic in the original ones. Such dual states were first 
introduced in Ref.-^, and were later defined differently in the three-qubit context by Borsten 
et.al.-ifi. These dual states have also made their debut to the physics of black holes admitting 
Freudenthal or Jordan duals^^ii^. 

Using these definitions we define the quadratic and sextic invariants as 



I, = -{A-D-B-C), 



(a ■ d — b ■ c) 



;i89) 



(The labelling convention and normalization for our invariants will be clarified below.) This 
form of the sextic invariant is deceptively simple. Its explicit form in terms of the dot 
product of Eq. fllSSp is 
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(190) 



We also recall that the explicit form of Ji is hiding its permutation invariance. A permutation 
invariant form is the one we used in Eq.( lll4p . i.e. we have ds^ = — /i(|A)) where |A) is the 
special state of Eq.f ll06p . Moreover, though the expression of J3 of Eq.f ll89p is similar to the 
one of Ji the invariant I3 is not invariant under the permutation of the qubits. 

Now we turn to the structure of quartic invariants. We have two of them and the simplest 
is the obvious expression 

h = DetA (191) 



i.e. the determinant of the 4x4 matrix of Eq.f ll84p . In order to present the definition of 
our last invariant we define separable bivectors of the form 



n 



^fia-^u(5 — A^^A 



a, (3, fi,^ = 1, 2, 3, 4. 



(192) 



Here our labelling convention A^q, indicates that a = 1,2,3,4 identifies the four-vector in 
question (i.e. A, B,C or D of Eq.f ll84p ). and the label /x = 1, 2, 3, 4 refers to the component 
of the particular vector. Now our last invariant is the quartic combination 



(193) 
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For the explicit form of this invariant we introduce the ■ product of two separable bivectors 

as 

{AAB)-{CAD) = 2{A ■ C){B ■ D) - 2{A ■ D){B ■ C). (194) 
Then the explicit form is 



6 



{A A B) ■ {C A D) + {A AC) ■ {B A D) - ]^{A ADf -]^{BA Cf 



(195) 



Let us now present the reason for considering these particular combinations for the set 
of algebraically independent S'L(2,C)^^ invariants. Let us consider the matrix 

n = A^gAg. (196) 

Then the characteristic polynomial of this 4x4 matrix is 

Det(lt -n)=t^- 4ht^ + Qht^ - Aht + II. (197) 

Clearly we have 

h = -Tin, h = ^[{Tinf (i98) 

J3 = l-[(jrnf - 3TrOTrfi2 + 2Trl]=^], (14)^ = Detl]. (199) 

This form of writing our invariants is related to the fact that there is a 1 — 1 correspondence 
between the 5*^(2, C) """^ orbits of four-qubit states and the 50(4, C) x 50(4, C) ones of 4 x 4 
matrices. 

The polynomial of Eq. fll97p first appeared in Ref .-41 its role as a characteristic polynomial 
has been emphasized in Ref.^^. The discriminant of this fourth order polynomial is the 
hyperdeterminanlr^ of the 2x2x2x2 hypercube Ajgjjj^jg. It is a polynomial of degree 
24 in the 16 amplitudes and has 2894276 terms^. It can be shown^Iiii that can be 
expressed in terms of our fundamental invariants in the form 

256L>4 = S^ - 27T^ (200) 

where 

S = {l! - l!) + - /1/3), T = {II - - h) + (Is - hhf. (201) 

In closing this appendix we briefly discuss some results on the full classification of en- 
tanglement classes for four qubits^*^. By entanglement classes we mean orbits under 
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SL{2,C)^^ ■ Sym4 where Sym4 is the symmetric group on four symbols. The basic result 
states that four qubits can be entangled in nine different ways^i^. It is to be contrasted 
with the two entanglement classes^ obtained for three qubits. 
Let us consider the matrix 

= (202) 

If A is the special matrix of Eq. fll06p used in the black hole context TZ/^ is just 2Vl of 
Eq.f llOSp . If the matrix TZ/^ is diagonalizable under the action 

f 83^82 \ 

7^Ah^57^A5-^ 5= , 5„g5L(2,C) (203) 

y Si ® So J 

we say that the corresponding four-qubit state |A) is semisimple. If TZ\ is nilpotent then we 
call the corresponding state |A) nilpotent too. It is known that a nilpotent orbit is conical 
i.e. if |A) is an element of the orbit then t|A) is also an element for all nonzero complex 
numbers t. Hence a nilpotent orbit is also a GL{2, C)*^^ orbit. (Recall that our V'^, is in 
the GL{2,C)^'^ orbit of the original V of Eq. fl62[ ) It is clear that for nilpotent states all 
of our algebraically independent invariants are zero. These are the states we associated to 
extremal black hole solutions of BPS and non-BPS type. 

A generic semisimple state of four qubits can always be transformed to the forn>^ 

\Gabcd) = ^^(10000) + 11111)) + ^^(10011) + 11100)) 

+ ^(10101) + 11010)) + ^(10110) + 11001)), (204) 

where a, b, c, d are complex numbers. This class corresponds to the so called GHZ class found 
in the three-qubit case^. For this state the reduced density matrices obtained by tracing 
out all but one of the qubits are proportional to the identity. This is the state with maximal 
four-partite entanglement. 

Another interesting property of this state is that it does not contain true three-partite 
entanglement. A straightforward calculation shows that the values of our invariants 
(Ji, /2, /s, 14) occurring for the state \Gabcd) representing the generic class are 



= i[a2 + 52 ^ ^2 ^ ^2]^ ^ 1 [(^^)2 ^ ^^^Y ^ (^^)2 ^ (^^)2 ^ (^^)2 ^ ^^^)2]^ (205) 
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= -[{abcf + {abdf + (acdf + (bcdf], h = abed, (206) 

hence the values of the invariants (4/i, 6/2, 4/3, J|) are given in terms of the elementary 
symmetric polynomials in the variables (xi, X2, X3, X4) = {a?, 6^, c^, d"^). On the generic class 
\Gahcd) the value of can be expressed as^^-i^ 

A = ^n,<,(x, - x,f = ^V{a\ b\ c\ d^, (207) 

where {xi,X2,X3,X4^) = {a'^,b'^,c^,d'^) and V is the Vandermonde determinant. 
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